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Weak turbulence of magnetohydrodynamic (MHD) waves in strongly magnetized plasma is stud- 
ied when the plasma pressure is less than the magnetic field pressure. In this situation the main 
nonlinear mechanism is the resonance scattering of fast magneto-acoustic and Alfvenic waves on 
slow magneto-acoustic waves. As a result, the former waves serve as the high-frequency waves with 
respect to the latter ones so that the total number of HF waves - an adiabatic invariant - conserves 
additionally. In the weak turbulence regime this invariant is shown to generate the Kolmogorov 
type spectrum with a constant flux of HF waves towards large-scale region. In the short-wave re- 
gion another Kolmogorov spectrum can be realized with a constant energy flux. The explicit angle 
dependences for both types of turbulent spectra are found for the propagation angles close to the 
direction of a mean magnetic field. 

I. INTRODUCTION 

The central place in theory of turbulence occupies the notion of spectrum of turbulence - the energy distribution 
on scales. The problem of its finding is one of the most difficult and nowdays is far from the complete solution. 
For developed hydrodynamic turbulence the pioneering works of A.N. Kolmogorov and A.M.Obukhov § of 1941 
about self-similar nature of turbulent spectra determined development of theory of turbulence for many years. In the 
seventieth years mainly by the efforts of V.E.Zakharov the Kolmogorov-Obukhov ideas were fruitfully applied to the 
weak wave turbulence theory (for details, see monograph || and also the first papers devoted to this subject j§-§[)- 
Wave turbulence in some sense occurred to be simpler than developed hydrodynamic turbulence. The main reason 
of such simplicity is connected with the wave dispersion when there exists such a region of wave intensities when 
the nonlinear interaction between waves can be considered weak in comparison with dispersive effects. If initially 
phases of waves are distributed randomly then the nonlinear interaction can provide a weak correlation in phases of 
the interacting waves. By this reason such ensemble of waves can be described in terms of pare correlation functions 
the Fourier spectra of which coincide, up to the multiplier, with a number of waves rik (occupation number) with 
the definite wave vector k. In their turn, the occupation numbers nk obey the kinetic equations for waves. The 
Kolmogorov spectra in this theory arise as stationary scale-invariant solutions of the kinetic equations annulating 
their collision terms. These spectra, unlike the thermodynamic equilibrium ones, can be related to the solutions 
with nonzero fluxes over scales. They realize a constant flux over scales of some integral of motion: energy, number 
of waves, etc., in the so-called inertial interval. It is important that if for developed hydrodynamic turbulence the 
notion of the inertial interval - the region where influence of both pumping and dissipation can be neglected - in fact 
represents the hypothesis of locality, for the weak turbulence theory this property can be checked explicitly. 

It is necessary to note that the main mass of works devoted to the Kolmogorov spectra of weak turbulence consider 
isotopic media (the corresponding bibliography can be found in H ) . Influence of anisotropy, for instance, magnetic field 
in plasma, has been studied in the less extent. The first such example of Kolmogorov type spectra was considered by 
the author in 1972 Q for weak turbulence of magnetized ion-acoustic waves in plasma. For this case the collision term 
of the kinetic equation was shown to be invariant with respect to two independent scalings along and perpendicular the 
magnetic field that allowed, by means of generalization of the Zakharov transformation, to construct the Kolmogorov 
spectra with the power dependences relative to longitudinal (k z ) and transverse (fc_i_) components of the wave vector. 
Later, the ideas of these papers were used for finding Kolmogorov spectra for drift waves in plasma and for the Rossby 
waves (see, for instance, PJTo|). 

The present paper is devoted to weak turbulence of magneto-hydrodynamic waves in strongly magnetized plasma 
when the thermal pressure of plasma nT is small compared with the magnetic field pressure H 2 /(8tt) (j3 — RnnT/H 2 <C 
1). In this case, in comparison with Q, turbulent spectra are defined from solution of three linked kinetic equations 
for Alfvenic, fast and slow magneto-acoustic waves. 

For f3 <C 1 the main nonlinear mechanism of the MHD waves is a scattering of fast magneto-acoustic and Alfvenic 
waves on slow magneto-acoustic waves (Section 2). In these processes, which can be considered as a partial case of 
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decays of one wave to two another waves and the oppsite process of fusion, Alfvenic and fast magneto-acoustic waves 
play a role of high-frequency waves against slow magneto-acoustic waves. In each act of scattering the frequency change 
of the former waves (we will call them A-waves) is small enough, due to smallness of the parameter (3. By this reason 
this process is familiar to the Mandelstamm-Brillouin scattering - scattering of electromagnetic waves on acoustic 
phonons. According to separation in time scales - devision of all waves onto HF and LF waves - the decay interaction 
conserves, besides the energy, the adiabatic invariant - the total number of HF waves. However, this analogy with the 
Mandelstamm-Brillouin scattering has not been exhausted by the said above. It turns out that the matrix element of 
this interaction has maximum for the maximal value of longitudinal projection of momentum transfered by A-waves 
to slow magneto-acoustic waves. In particular, this result can be extracted from the expression for the growth rate 
of decay instability for the monochromatic Alfvenic wave obtained in 1962 by A.A.Galeev and V.N.Oraevskii p2[ . 
Remind that for the Mandelstamm-Brillouin scattering the matrix element is also proportional to square root from the 
value of transferred momentum, that provides maximum for the back scattering of electromagnetic waves. Because 
of such behavior of the scattering amplitude of A-waves it is natural to assume that a stationary distribution of 
waves over angles will be very anisotropic, concentrated in the fc-space along the magnetic field direction. Under such 
assumptions the kinetic equations possess two additional symmetries - invariance with respect to two independent 
stretching along and in transverse direction to the mean magnetic field that allow one to use the transformations of 
the paper . Due to two these symmetries of the kinetic equations it turns out that in the transparency region there 
are possible two scale-invariant (against longitudinal and transverse wave vectors) Kolmogorov spectra corresponding 
to a constant flux of energy to the short-wave region - the direct cascade and to a constant flux of A-waves towards 
the large-scale region - the inverse cascade. This paper is based on the old results of the author |l3| published in the 
form of preprint in Russian and unknown by this reason abroad (that is natural) . But it turns out that these results 
are also unknown for Russian readers. In spite of more than twenty five years history nobody has not repeated these 
results. It is necessary to note, however, that recently the question about MHD turbulence has been considered in 
another limit (5 3> 1 in JT3]. This limit differs significantly from that considered in the present paper. First, at j3 ^> 1 
plasma can be considered as almost incompressible fluid and, secondly, in this limit there is no essential difference 
between Alfvenic and slow magneto-acoustic waves: the latter waves have the same dispersion law as Alfvenic waves 
differing from them by polarization only. Such degeneracy sufficiently changes the character of nonlinear interaction. 
In spite of these facts in this situation two kinds of Kolmogorov spectra are possible with the same dependences on 
wave numbers as the obtained ones in the present paper. However, the physical motivation of existance of the two 
Kolmogorov spectra at (3 » 1 is different. 

Plan of the paper is the following. In Section 2 we introduce canonical description of an ideal MHD following to the 
original paper of 1972 by V.E.Zakharov and the author jll]] and their recent review |l2] ]. By means of the Hamiltonian 
description in the next - third section the average equations are derived for A-waves with account of interaction with 
slow magneto- acoustic waves. It is shown that from the side of A-waves to slow plasma motion there acts the HF 
force. The potential of this force is negative. Therefore, unlike the interaction of Langmuir waves with ion-acoustic 
waves Eqj, plasma is drawn into regions of A- wave localization forming there density humps. In the same section 
we analyze stability for the monochromatic A-waves. Section 4 is devoted to the Kolmogorov spectra of weak MHD 
turbulence. 



II. VARIATIONAL PRINCIPLE AND CANONICAL DESCRIPTION 

Consider the equations of ideal MHD for barotropic flows of plasma when the internal energy of plasma e depends 
only on its density p: 

g+divpv-0; (2.1) 

^ + (vV)v = -Vw + -!- [rotH x H] ; (2.2) 
at Anp 

U — = rot [v x H] . (2.3) 

Here v is plasma velocity, w entalphy connected with pressure p = p{p) and internal energy e by the relations: 

j dp d 
dw = — ,w= -^-e{p). 
p dp 
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Let us formulate the variational principle for this system. 

First of all, it should be noticed that as follows from the equations ( O 2.3 ) the vector H/p is advected by fluid 



particles, by another words, each magnetic line moves together with its own particles. This is the well known fact 
of frozenness of magnetic field into plasma (see, e.g. j25|). The given circumstance allows one to consider both the 
magnetic field H and the density p playing the roles of generalized coordinate. 

To formulate the variational principle we shall use the well known expression of the Lagrangian of electromagnetic 
field and particles (fluid) within pl| . The Lagrangian L should be written in the MHD approximation. In particular, 
this means that one needs to neglect by the contribution from the electric field E in L in comparison with that from 



the magnetic field since E ~ v/c H <C 1 (here c the light speed). Secondly, we shall account Eqs. (|2.l|), ( f2.3| ) and 
div H = as constraints. As a result, the Lagrangian can be written as follows 

£?!- e ^-5! + S.rS- rot rvxHiWsfS 



L = —- e (p) - g^T + S • {-qT - rot I v x H]J + + div pvj + V div H. 

Here S, $ and ip are unknown Lagrange multipliers depending on r and t. 
Next, using this expression of the L we write down the functional of action: 

I = J Ldt dr, 

which variations relative to the functions v, p and H yield the equations: 

pv = [H x rotS] + pV$ , (2.4) 



~dt 



(W)$ 



+ w(p) = , 



(2.5) 



-TTT + ~ A [v x rot S] + = . 

at An 



(2.6) 



The first equation in this system gives the change of variables: velocity v is expressed in terms of new variables S 
and $. It is necessary to note that this change of variables is ambiguous: it is possible to add S the vector So, and 
to $ the scalar $ satisfying the equation 

[H x rot S ] + pV$ = 0. 



Two other equations in this system - Eqs. (2.5) and ( |2.6| ) - represent themselves the Bernoulli equation for the 
potential and the equation of motion for a new vector S which contains unknown potential ip. The potential -0, in 
turn, is defined by fixing gauge of the vector S. For example, for the Coulomb gauge (div S = 0), ip is determined up 
to arbitrary solution ipg of the Laplace equation Aipo = 0: 

tp = — div [v x rot S] + %/jq . 

In particular, if v — > 0, H — * Ho, p — > pa at infinity on r then the value tpa is convenient to be chosen so that S 
at r — * oo. Then obviously ipo = — (Hq • r)/(47r). 







Now one needs to chec k th at the system of equations (2.4-2.6) does not contradict to the original system of MHD 
equations. By plugging (2.4) into the equation of motion (2.2), after simple transformations, it is possible to verify 
that Eq. ([2.2]) can be written in the form: 



V dt 



.9$ v 2 

( vV )$- — + w (p) 



H OS 

7 X rot a 



H 

47T 



[v x rot S] 



0. 



According to (2.5,2.6) the obtained equation satisfies identically. Thus, one can say that the new system 
(2.1 2.2,2. 5, E^q) is equivalent the MHD system. Really, due to the formula (2.4) any solution of the system 
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( f2~l 2.5 ,2.5 JT6]) generates some solution of the original MHD system. Under assumption about uniqueness of the 
Cauchy problem for the systems (2.1-2.4) and (2.1,2.3,2.5,2.6) opposite statement is valid also: for any solution of the 
system (2.1-^4) is possible to put in the correspondence some class of solutions of the system (2.1,2.S.2.5,2.6). For 
this case it is enough by using a set of quantities v, H, p at some moment of time to one to construct various sets of 
quantities S and <£, satisfying the formula ( [2.4| ) and to take them as initial conditions for the system (2.1,2.3 2.5,2.6). 

After that, by means of the Lagrange function, we determine the generalized momentum and construct the Hamil- 
tonian of the system, by the standard way: 



H = J ((S • H t ) + <&p t — L 



- L)dr = 



■1 Jj2 

— + £ (P) + — - ^ div H \dr, 



which coincides with the total energy of the system. The equations of motion (2.1,2.3,2.5 2.6) is this case are nothing 
more than the Hamilton equations: 



dp _SH a$ SH 



(2.7) 



<9H 

~dt 



SH 
~SS' 



as 

at 



SH 
SB. 



Respectively, the variables (p, <E>) and (H, S) are pairs of canonically conjugated values. Change of variables (2.4) and 
canonical description ( [2.7|) were introduced for ideal MHD by Zakharov and the author in 1970 in the paper |ll| . 

The transformation (2.4) represents analog of the Clebsch representation in the ideal hydrodynamics. Respectively, 
the fields H and S in the formula ( |2.4D play the same role as that of the Clebsch variables (about Clebsch variables see 
p8[ , [ p6| , and recent review |l^]). Later, in 1982, the authors of Ref Jl7|]came, in fact, to the same change: the velocity 
vector and magnetic field were represented through the scalar Clebsch variables that after simple transformations can 
be reduced to fl2.4p. 



MHD flows describing by means of (2.4) as well as flows of ideal fluids parameterized by the Clebsch variables are 
related to the partial type of flows. For example, for such MHD flows the topological invariant of linkage of magnetic 
field and vorticity lines, the so-called cross helicity |l(| / — J(v • H)dr is identically equal to zero. 

In 1995 Vladimirov and Moffatt [15[ for ideal MHD found the analog of the Weber transformations: 



v = it fe(a)Va fe +V$l-[Hx rotSl. 

P 



(2.8) 

r(a, t) which defines particle 



Here a = a(r, t) are Lagrangian markers of fluid particles (this mapping is inverse to r 
trajectory with the marker a ), u (a) is a new Lagrange invariant. 

The Weber transformation (2J3) is the general transform which, in particular, contains at Uo = 0, the change (2.4) 
to what the authors of [fl5| did not pay attention. It is interesting to note that the equations of motion for potentials 
$ and S for the general transform ( |2.8| ) have the same form as for (2J5) and ( 2.6). I f $ and S are equal to zero at 
t = 0, then Uo(a) is nothing more than the initial velocity. Just the first term in (2^) provides nonzero value of the 
topological invariant /. One should note also that this term, as for ideal hydrodynamics, is nonlinear if one proceeds 
expansion over small amplitudes. Recently, Ruban [[l9| (see also [^0|) clarified the physical meaning of the vector S. 
Curl of the vector S can be expressed through the shift d between electron and ion (as fluid particles) in the point r 
at the moment of time t if initially their coordinates coincide: 



rotS 



e , p(r,t) 
Mc p (a) ' 



Here M, e are ion mass and its charge, respectively, Po(sl) the initial distribution of plasma density. 

Introducing the canonical variables allow by the standard way (by means of the perturbation theory against small 
amplitudes of waves) both to classify and investigate all nonlinear processes. To this aim one needs in the expression 
for the velocity ( |2.8[ ) as well as in the internal energy to perform expansion in powers of the canonical variables. In 
the presence of the external homogeneous magnetic field Ho, in the linear approximation, one needs to keep the linear 
terms relative to $ and S, neglecting by the first (nonlinear) term in (2.8). As the result, in the velocity expansion 



v + Vi 



(2.9) 
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the first order term is written as 



v = — [H x rot S] + V$. 

Po 

The three independent pairs (div H = div S = 0) of canonically conjugated quantities will correspond to three 
types of waves. In the linear approximation, obviously, these waves will not interact. Their dispersion laws and 
polarizations can be found from analysis of the quadratic (relative to the canonical variables) Hamiltonian Ti.Q. 
Three- wave interaction will correspond to the cubic (relative to the canonical variables) term TL 3 . Its value will be 
defined from the quadratic (against wave amplitudes) additions to the velocity 

vi = ^i[H x rot S] + — [h x rot SI. 

Po Po 



In this expression we take into account only 'wave' degrees of freedom and neglect by the first term in (2.8). Here h 
and pi are the magnetic field H and density p variations from their equilibrium values Ho and po, respectively. As a 
result, the Hamiltonian represents a series with respect to powers of the wave amplitudes: 

H = H +H 3 + (2.10) 

where the quadratic Hamiltonian is 

' Po^l h 2 2 p\ 



and 



n 3 = J {po(v • v x ) + ^-vl + qc 2 s ^dr 

2 



is the cubic Hamiltonian. In these formulas square of the sound speed c s and the dimensionless coefficient q appeared 
from expansion of the internal energy in powers of pi : 
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Let us now perform the Fourier transform for coordinates and then introduce new variables dj{k) (j — 1,2,3) by 
means of the following formulas 

h(fc) = ei(fc)V^i(ai(A;) + a^(-fc))+ (2.11) 



+e 2 (fc) V A z V27rw/(a/(fc) +a*(-k)); 



1=2,3 



S(k) = - ie i(fc)-_(ai(fc) - ol(-fc))- 



1=2,3 



P^= E fe) ^M) + a*_ k {l)) 

1=2,3 \ C S / 



1=2,3 



\ 2po^; 



1/2 



Hi(ai(k) - af(-fc)). 
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Here 



w x {k) = |(k- V A )|, 



k 2 V\ + k 2 c 2 + 2(k • V A )kc s ± \k 2 V 2 + k 2 c 2 s - 2(k • V A )k 



are the dispersion laws of Alfvenic waves (1), fast (2) and slow (3) magneto- acoustic waves, respectively; 



ei(fe) 



[k x n ](k ■ n ) 
[k x n ]||(k • n ) 



ea(fe) = 



[k x [k x n ]] 
fe|[k x no] | 



are corresponding to them the unit polarization vectors; no = Ho /-Ho the unit vector along mean magnetic field; 
= Ho/(47rp ) 1/2 the Alfven velocity; 



-A*3 



uj-k 2 c 2 s 
lu 2 — k 2 c 2 



1/2 



-1/2 



The given change to the new variables Ofc(j) represents the canonical U — V transformation diagonalizing the Hamil- 
tonian TIq: 



Ho = J2j u 3 {k)a {k)a*{k)d\<L. 



In this case the amplitudes cij(k) have the meaning of normal variables, and, respectively, the equation of motion 
have the standard canonical form: 



da {k) 
dt 



. SH 



In the linear approximation Oj (fc) obey the equations: 



daj(fc) 
dt 



+ iujj{k)aj(k) = , 



namely, with time the amplitude modulus |aj(fc)| does not change, but the phase increases linearly with t. 

In order to find the exp ressi on for the Hamiltonian of interaction in terms of the a,j(k) variables one needs to 
substitute transformation ( [2.11 ) into ( 2.10 ). As a result, the Hamiltonian of interaction will represent the integro- 
power series relative to theses variables. In the lowest order with respect to the wave amplitudes the main nonlinear 
process will be the resonant three-wave interaction. The corresponding Hamiltonian is equal to 



«w = \J J2[ V ^ a l( k ) a m(h)a n (k 2 ) 



Z2 dkdkidk2- 



(2.12) 



This Hamiltonian can be obtained after substitution of the transform (2.11) into the cubic Hamiltonian Ti^ and 
subsequent extraction from there the resonant terms. The rest terms in H3 are small: they can be excluded by means 
of the canonical transformation (for details see the review ]1^] ) . One should note that calculation of matrix elements 
in this scheme is a pure algebraic procedure requiring performance of the Fourier transform in all integrals, 



kki hi 



substitution (2.11) and forthcoming symmetrization with respect to the ak(i) variables, for instance, in (2.12) against 
pairs (ki,m) and (k2,n). 
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III. AVERAGE EQUATIONS 



Expressions for dispersion laws and matrix elements of interaction can be sufficiently simplified for plasma with 
low j3 — %nnT / H 2 (it is ratio of thermal plasma pressure nT and magnetic field pressure H 2 /8tt. The condition 
/3 <C 1 means that Va 3> c s . In this limit fast magneto-acoustic waves have isotropic dispersion law u>2 = IzVa and 
their phase (as well as group) velocity coincides with the value of the group velocity of Alfvenic waves. In the linear 
approximation the plasma velocity in Alfvenic and fast magneto- acoustic waves is given by the expression 

v HF = —[Ho x rot S]. 

Po 

The potential part of the velocity V<3? occurs to be small due to smallness of the parameter (3. For slow magneto- 
acoustic waves, on the contrary, the main contribution is given by the potential part - the velocity turns out to be 
directed along the magnetic field H : 

v s = n ^-, (3.1) 

oz 

and the dispersion law for slow magneto-acoustic waves becomes strongly anisotropic: 

uj 3 = Q a = \k z \c s . (3.2) 

The transverse components of the velocity for these waves [Ho x rot S]/po are compensated by Vj_<!>. 

If plasma is collisionless and strongly isothermal (T e 3> T,), the slow magneto- acoustic waves represent themselves 
magnetized ion-acoustic waves (for details see ||). In this case in ([^) c s — ^/T e /M. 

As far as the nonlinear interaction of the MHD waves concerns, for strongly magnetized plasma the main nonlinear 
process is the process of scattering of the Alfvenic and fast magneto-acoustic waves on the slow magneto-acoustic 



waves (that can be verified directly by comparing the computing matrix elements V lmn for (2.12)). In this process 
the former waves (further we shall call them as A-waves) plays the role of the high-frequency (HF) waves relative 
to the latter ones (these waves shall be simply called sound or S- waves). This conjecture follows directly from the 
analysis of the resonant conditions for the given type of decay: 

uo A {k) = u A (ki) + n»(*a), k = ki + k 2 . (3.3) 

Qualitatively it is easily to understand how this interaction looks like. While propagation of the packet of A-waves 
the mean characteristics of plasma (its density and mean velocity) due to the action of A-waves will be slowly varied. 
By this reason the mean Alfven velocity will differ from its local value by the quantity AVa = — VaPis/^Po) where 
Pi s is low- frequency (LF) density variation. It results in the frequency addition of A-waves Alo p ~ fcAV^. Due to 
slow motion of plasma with the drift velocity Vjj the frequency of A-waves changes at Aw^ ~ kv£, . The ratio of these 
two additions, Awp and Aw p , however, occurs to be small: ~ c s /Va- Thus, the main interaction is the scattering on 
the LF density fluctuations. At the same time the LF plasma characteristics will be changed due to the action of the 
HF force induced by A-waves. 

The most simple way to find the expression of the HF force is to perform average of the Hamiltonian over the HF 
oscillations. The result of this average is the following 

n = n + n iat , (3.4) 

where 

'W{>**-*">+g}* + /{ a £ + 4&}*, 




Here angle brackets mean average over high frequency. The first integral in TIq corresponds to (linear) A-waves, the 
second one describes magnetized acoustic waves and the last term is responsible for nonlinear interaction. Variation 
of the Hamiltonian of interaction on p\ s gives the expression for the HF potential : 
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In according with (3.5) the equation of motion for the potential Q s takes the form: 

d<P s 2 Pls ([H x rot S] 2 ) 

"ar + c ^ = — M • (3 - 6) 



It is important to notice that the HF potential ( |3.5| ) is negative. This means that in the region of localization of 
A-waves the HF force will form, instead of density wells, as it is for the interaction between Langmuir and ion-acoustic 
waves (see 0), the density humps. 

The equations of motion are closed by the continuity equation for p\ s which, in accordance with (EOj), has the form: 



From (3.6) and (3.7) we have 



dp ls <9 2 $ s 



d 2 Pis 2 9 2 Pis 1 d 2 2 

-dfi ~ C ^ = -2p- d^ {[Ho X r0t S] >' (3 - 8) 



To write the equation for A-waves one needs to make average explicitly in the Hamiltonian of interaction TLint- It 
corresponds to keeping in 7ii n t terms containing products a A a Al where the index A = 1, 2 enumerate the HF waves: 



Here 



flint = - ( Pl o R ^F^al(k)a Xl {ki)S k - k ^ K dkdkdK. 

= K(fc)c Al (fc 1 )) 1 / 2 (n A (fc) • n Al (fcx)), n 2 = ^, m = -[n 2 n ] 

As the result the equations for A-waves have the form: 

da x (k) . SH mt 

gt + iu)x{k)a\{k) = -i j-^-j ; A = 1,2. (3.9) 

For isothermal collisionless plasma (T e w Tj) slow magneto- acoustic waves are absent due to strong Landau damping 
on ions. Correspondingly, the decay interaction of A-waves transforms into the induced scattering on ions. In this case 
the equations ( |3.8| ) have to be changed by the drift kinetic equation for slow variation of distribution functions 
of ions fi (compare with [p9[): 

dfi , df z 1 d df 
ot oz M oz av z 

together with the quasi-neutrality condition for slow motions (f^ = k z cs *C u) v i) 

8ni = y /»^v = Y ei P = ^ ( 3 - n ) 

where f is the equilibrium distribution function of ions, and (p the LF electrostatic potential. In this case the equation 
of motion for A-waves retrains the form of (|3.9[), and the density is expressed linearly through the HF potential by 
means of the Green function for the system ( |3. 10 3.11 ): 

_ pis(n,n) n K 2 e e e t 

Lr K f2 = Tj = 2 i ■ l 3 ' 1 ^) 

U K n lo^ e e + Ei 

Here Pi s (k,Q) and U K si are the Fourier images of the LF density and the HF potential, respectively, e e ,i partial 
dielectric constants of electrons and ions which are equal: 

1 



O 1 

d 
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47re 2 [K z (df /dv z ) 

-ctv . 



Mk 2 I Q — k z v z 



where r\ — T e / (Airnoe 2 ) is square of the Debye radius. 

In non- isothermal plasma (T e 3> Tj) the Green function ( |3.12 ) transforms into 



Lr K f 



O 2 



that coincides with the expre ssion g iven by the equation (3.8) 



The system of equations ( p.lC| )-( 3.12 ) completely describes interaction of A- waves in magnetized plasma with 
arbitrary ratio of ion and electron temperatures. In this case, however, the Hamiltonian Ho +-ffjnt ^ s n0 ^ conservative 
quality due to the Landau damping on ions. 

IV. INSTABILITY OF MONOCHROMATIC WAVE 

Let us now analyze the obtained equations. We start from study of dynamics of the narrow packet of A- waves. 

A qualitative understanding about this process can be obtained from stability analysis of monochromatic A-wave. 
In the sake of simplicity we shall restrict by consideration of stability of Alfvenic wave in the hydrodynamic limit. 
For collisionless plasma the latter assumes that the phase velocity of bending Q/k z for A- waves exceeds the thermal 
ion velocity VTi- In this case for slow motion one can neglect by the Landau damping on ions and use the equations 



(3.8) or (3.12). One should remember that in strongly non-isothermal plasma plasma the sound waves are the eigen 
oscillations but at the same time in plasma with T e w Ti sound waves represent the induced oscillations of the plasma 
density. However, at CI/k z 3> vti the hydrodynamic description can be applied in both cases. 



Then it is convenient to express p\ s through the normal variables aa(k) = b k : 



The equations of motion for b{k) are obtained from variation of of the full Hamiltonian Hq + H^ n + 

db k 



;jt +iSl(k)b k = -iSH hA /Sbt. (4.1) 



In the equations (|3.9| ), flLl| ) the solution 

A 

a\(k) = — ^Sxie~ luot 5 k -k ^ b k = 0, lu = wi(fc ) 

corresponds to the monochromatic Alfvenic wave. The amplitude of Alfvenic wave here is chosen by such a way so 
that | A | 2 coincides with the energy density of oscillations W. 

Linearizing of the equations ( 3.1C| )-( [3.12 ) on the background of exact solution and assuming for perturbations: 



Sa x {k) ~ e - l(n+ " o)t 4- fco - K , 

Sa* x (k) ~ e-^-^t^+K, 
for Q we have the following dispersion relation: 

WG | l^q-o+J 2 ! \ F koM-^ 2 1=1 (4 2) 

We shall present now results of investigations of the dispersion equation ( |L2| ) in the different cases in dependence 
on the energy wave density W and on the temperature ratio. 

At T e 3> Ti and sufficiently small amplitudes the decay instability takes place with generation of ion magnetized 
sound 122] . For this instability the eigen frequency f2 is expressed through the matrix element of the decay interaction 
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y k fx = [^) 1/2 ^ (4.3) 



and the quantity W: 



1. » . fir ,. , ~, M 2 1A 



1/2 

n = J [wo - wa(*o -k) + fi(«)] ± | J [w - w A (fc„ - «) - fi(K)] a - ^|^ fco _ K , K | 2 | • (4.4) 

Hence it follows that the instability takes place near the resonant surface 

uj = uj\(k — k) + fi(«) (4.5) 



with maximum of the growth rate 



w n Hlx 



8nT w '~ k ° k ° 



1/2 



(4.6) 



The growth rate width as a function of frequency occurs to be of the order of magnitude of the maximal growth rate 



(4.6) 



Because the matrix element is proportional to the frequency of slow sound the maximum value of the growth rate 
on the resonant surface ( |4.5| ) is attained at the maximal value of |k z |. For decay on Alfvenic wave and slow sound 
max |ac z | Ri 2|A:o z |, so that the secondary Alfvenic wave propagates in the opposite direction to the pumping Alfvenic 
wave. Such behavior of the decay instability is typical for the Mandelstamm-Briullien scattering, the matrix element 
of which is proportional to square root from the sound momentum transfered by scattering light. For light such 
dependence provides the maximal back scattering. 

It is not difficult to investigate the decay instability for all o ther possible channels of decay A — > A + S. The growth 
rate in all these cases are of the same order of magnitude as ( |4.6| ) : 

T - (Lu Q s W/nT)^ 2 . 

This instability takes place at 

W/nT < P 1/2 . 



With increase of W the decay instability is transformed. At W/nT > /3 1 / 2 in the dispersion relation (4.2) it is 
possible to neglect by £l 2 against fl 2 . Then the unstable wave vectors will lie on the surface uJi(ko) = uJ\(ko — k). This 
instability is called as the modified decay instability [p8|]3C| l . For the interaction of Alfvenic waves and slow sound this 
instability has the growth rate maximal at k z = 2ko z : 

V3 ( W \V3 
-w 



fa* ' (47) 



2 "VoVi 

Value of this growth rate does not depend on temperature and therefore this instability takes place at W/nT > 1 up 
to the values when the main approximation - adiabaticity approximation - looses its applicability: T ~ luq. 
For another channels the instability with growth of W/nT has the same character: at W/nT > /3 1 / 2 the growth 



rate is maximal in the re gion k ~ fco an< 4 has the same order of magnitudes as (4.7) 



The decay instability (4J3) for arbitrary ch ann el A — > A + S, as it is easily seen, relates to the convective type of 



instabilities. The excited waves, according to (4.4), have the group velocities strongly different from the group velocity 
of the pumping wave. Therefore for the wave packet with the characteristic scale L this instability will be essential 
only for large enough lengths L when the amplification coefficient G exceeds a value of the Coulomb logarithm A: 

G = TL/Va w A. 

For less lengths L, the decay instability is not important: during propagation of perturbation of through the whole 
packet perturbations amplify for a small value. In this case dynamics of the packet will be defined by slow processes. 
Among them the most important ones are such processes for which unstable perturbations propagate together with 
the wave packet. If it is a decay instability, then it has to be absolute (in the frame moving together with the packet). 
In particular, this is one of the reasons of appearance of collapse for fast magneto-acoustic waves and of affect of 
sound collapse on the fine structure of collisionless shocks in plasma |^3],^4| . Collapse of fast magneto-acoustic waves 
appears due to three- wave interaction in which only fast magneto-acoustic waves take part. 



10 



V. KOLMOGOROV SPECTRA 



In the previous section we considered the stability problem for narrow in fc-space wave packet. In this case for decay 
of monochromatic wave under the resonant conditions ( |4.5| ) (namely, for the maximal value of the growth rate ( |4.6] )) 
sum of the phases of exciting waves 4>a and <p s are strongly connected with the phase of the pumping wave 4*$: 

4>0 + 7r/2 = <t>A + 4s- 



(It is easily to check that this phase correlation is lost with leaving the resonance (4.5).) Simultaneously, a difference 
in phases for the pair of exciting waves with fixed k remains arbitrary. Both these factors introduce to the system of 
interacted triads, connected with the pumping wave an element of randomness. Thus, each triad is characterized by 
one random phase. At the next step - at the second cascade a new random phases are added so that a memory about 
the pumping wave will be lost. For multiple repetition of of this process the system of waves must transform in the 
turbulent state when the phases of waves can be considered random. Therefore the stochastization time should be a 



few inverse growth rate (4.6) 



Such scenarion of transition to turbulence seems to be sufficiently plausible. Now there are being performed a 



number of numerical experiments to check this hypothesis (see, for instance, ]33|], |32|). 

Hence it becomes clear that the regime of developed wave turbulence should be characterized by the wide spectrum 
of waves. For small enough intensities of waves it is enough one to restrict by consideration of the pair correlation 
functions only which obey the kinetic equations for waves. This regime is called weak turbulent. 

In the case of the weak MHD turbulence at (3 <C 1 we have three pair correlation functions defined by the following 
formulas: 

(a A (fc)a^(fci)) = N^Sx^Sk-k!, {b k b* ki ) = n k S k - kl 

where the quantities N£, n k , having a meaning of the occupation numbers, satisfy the following system of kinetic 
equations: 

n k = 2n J \V kl k 2 k\ 2 (N kl N k2 - n k N kl + n k N k2 )5 k+ k 1 -k 2 5n+u 1 -u 2 dkidk 2 , (5.1) 
N k = 271-/ \V kklk2 \ 2 (N kl n k2 - N k n k2 - N k N k , -fe<L-u;, -0,^1*2 (5.2) 



-2-7T / \V klkk2 \ 2 (N k n k2 - N kl n k2 - N k N kl )S kl - k - k2 S UJl - u -n 2 dkidk2. 

Here u> = u>(k), u>\ = u)(ki), and so on. In these equations (as well as below) we omit summation over A. In order to 
include it one needs to change dk\ — ► dk\, N k — ► N£, u k — ► tu k x, V kklk2 — > V k k \ an< ^ so on - 



The equations (5.1,^2) assume weakness of the nonlinear interaction between waves. In this concrete case the most 



essential criterion is 



n s > l/r, 



where r is the characteristic nonlinear time defined by the kinetic equations (5.1,5.2). To estimate the value of r one 
needs to take into account that in each act of decay and inverse process - merging of waves, the frequencies of A- waves 
change for the small value Auja = <C loa, namely, the energy transfer of A- waves along spectrum has a diffusive 
character. Due to this fact, we have the following estimate for r: 

1 W 

Notice, that this value for r exceeds significantly the stochastization time, defined by the inverse growth rate ( |4.6| ) 

r- 1 . 

Hence we have finally the criterion: 

w oX/2 



2 
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Next, let us include into the kinetic equations ( |5.l| ), ( |5.2| ) both sources of turbulence and its dissipation. For this 
aim in the left hand sides of the equations we introduce new terms T k Uk and "f k \Nk\, respectively. We suppose here 
that the pumping region (r k ,jk\ > 0) and dissipation region (1^,7^ < 0) are well separated in fc-space by the 
intermediate region - the inertial interval, where dynamics of turbulence is defined by nonlinear interaction between 
waves only. In the inertial interval we shall neglect by influence of both pumping and dissipation (that is necessary 
to be proved) then distributions rik and N k \ do not depend on the concrete form of 7^ and IV 

We would like to remind that in the theory of developed hydrodynamic turbulence to determine turbulent spectrum 
- distribution of energy for velocity fluctuations it is enough to use two hypothesizes of A.N.Kolmogorov 0. The first 
hypothesis about self-similarity says that the spectrum of turbulence is defined by the unique quantity P - a constant 
flux of energy along scales (from large to small ones where dissipation due to viscosity becomes essential) . The second 
hypothesis assumes that the interaction of fluctuations with different scales has a local character. 

If one applies these hypothesizes to the given case then spectra of turbulence in the in ertia l interval can be found 
by means of the dimensional analysis. In the given situation the kinetic equations (5.1), (5.2) have two conservation 
laws: conservation of the total energy and the total number of HF waves. To each integral of motion there should 
correspond its proper Kolmogorov spectrum. So, for a constant flux of the number of HF waves N£ 



we have the spectrum: 



For a constant flux of energy 



N kx ~ f#V 4 , n k ~ P]i 2 k-\ 



d { 



(5.3) 



one can get the estimate: 



N, 



pl/2 fc -3/2 



pl/2 fe -3/2 



(5.4) 



From conservation in the inertial range of the total number of HF waves and the energy it is easily to establish that 
the flux of HF particles is directed towards small k region, and energy flux is directed to the short wave region. 

These - sufficiently rough - estimations for the spectra (5.3), (5.4) can pretend only to the right dependence on both 
wave numbers and fluxes, but they don't account the diffusive character of decays. Notice also that these estimates 
are significantly based on an assumption of the interaction locality. 

The spectra ( |5.3| ) and (5.4) don't account also fine properties of distribution functions - their angle dependences, 
namely, they arc defined up to arbitrary functions of angles. To find these dependences one needs to solve the exact 
equations (5.1), fl5.2|) . It turns out that solution of these equations can be found for interaction of Alfvenic and slow 
acoustic waves (N2 = 0). For this case it is convenient to represent the equations (pj]), (5.2) in the form: 



U k2 \kk 1 Tk 2 \kk 1 dkidk2 , 



(5.5) 



N k = J (U k \ klk2 T k \ klk2 — Uk 1 \kk 2 Pk 1 \kk 2 )dkidk2 
where the following notations are introduced: 

Uk\k x k 2 = 27r|F fe , fe / fe J Sk-kx-kiSu-m-ws > 



(5.6) 



Tk\k x k 2 = N kl n k2 - N k n k2 ~ N k N kl . 
It is easily to see that the equations (pTSj) , (|5.6|) have the thermodynamic equilibrium solutions: 



N 

LOk +/x' 



T 
f2i- 
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- the Rayleigh- Jeans distributions which annulate the collision terms. 

To find the non-equilibrium distributions it should be noted that the function U has the following properties, (i) 
U is bi- homogeneous function of its arguments k z and k± with homogeneity degrees, respectively, equal to +1 for k z 
and —2 for k±. This means that if one performs stretching of k z , ki z , k2 Z in A times then U multiplies in A +1 times: 
U — ► \ +1 U. If one makes the similar transform for all k±: k± — > fikx, then U transforms as U — > fi~ 2 U . Besides, (ii) 
U is invariant with respect to rotation around z-axis - the direction of mean magnetic field Ho. 

Due to these properties solution is naturally sought in the form: 



n k = Ak%kP x , N k = Bk^y- A 



Consider the stationary equation (5.6) 



(Uk\k 1 k 2 Tk\k 1 k 2 ~ Uk 1 \kk 2 Tk 1 \kk 2 ) — 0- 



(5.7) 



(5.8) 



Let us make a mapping of the integration area of the second integral (which is given by its resonance conditions 
- (5-functions) to the integration area of the first integral. For this aim it is convenient to introduce the complex 
variables £ = k x + ik y . Then the integration area defined by the corresponding conservation laws 



0, 



Ci - C - C 2 = o, 



u>i — UJ — 0,2 = 0, 
with the help of the mapping relative to all k z and £: 



h, ry 



(5.9) 



Ki = k z ^, Ci = C^7, 



k z -? = k'i — , 



(,2 - C"^, 



transforms into the integration area of the first integral in (5.S). Each such transformation (separately with respect 
to k z and £) represents itself the operation of inversion: for z-components of wave vectors against the point k z , and 
for transverse components relative to the circle with radius \k±\. Under this mapping the vector k transforms into 
ki, ki into k and k2 into kj. Simultaneously all z-components are stretched in the \k z jk z \\ times, and all transverse 
components get the factor \k±/ki±\. Besides, the rotation on the angle arg(£/£i) around z-axis takes place. 
As the result, due to the properties of both U and T, 



U kl \kk 2 -> \k z /k zl \ +1 \k ± /k 1± \ 2 U k \ klk2 , T kl 



\kk-i 



\k z /k zl \ 2a \k ± /k 1± ^T k{klk2 , 



the integrand in (5.8) is factorized: 



Uk\k x k 2 Tk\k x k 2 



1 - 



/ k z \ 2 «+ 4 / fc_L \ 2/3+4 

\k~~J \~kTj 



dk\dk2 = 0. 



Hence it follows that, besides the thermodynamic equilibrium spectra (which annulates T), the following solution is 
possible: 

n k = Ak; 2 kf, N k = Bk~ 2 k- 2 . (5.10) 
These spectra correspond to the solution which was obtained previously from the dimensional analysis for the constant 



flux of HF waves Pn- Connection between the coefficients A and B in (5.1C) is determined from solution of the 
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stationary {d/dt — 0) equation ( |5.l| ). Hence one can get that for this case total energies containing in Alfvenic waves 
and in slow magneto- acoustic waves are of the same order of magnitude: c s A ~ VaB. 

It is worth to note that the whole set of the transformations ( |5.9D forms the group G. This group is direct product 
of two groups G(l) and G(2): G — G(l) x G(2). The group G(l) acts in one-dimensional space (£; z -space), and G(2) 
in the two-dimensional one (/c^-space). These transformations allow one to factorize the collision terms. First these 
transformations (in the ID case of the frequency space) were found by V.E.Zakharov |3]-[7|]. The generalizations of these 
transformation to the both 2D and 3D cases for isotropic models were introduced by A.V.Kats and V.M.Kontorovich 
in 1970 jn). The transformations (5.9) were found by the author in 1972 
so-called quasi-conformal transformations. 



They represent a partial type of the 



To find another non-equilibrium solution of (5.4) it is convenient to introduce the energy density in the fc-space 
Ek = UkNk + VL^nu- From (5.5) and (5.6) follows that this quantity obeys the equation: 



de k 
dt 



-TjT" — / \ u >kUk\k 1 k 2 Tk\k 1 k 2 ~ UJ k Uk 1 \kk 2 Tk l \kk 2 ~ ^kUk 2 \ kl kTl 



dki dk2 



(5-11) 



Consider stationary solution of this equation. As before, we shall seek for solution of ( |5.11 ) in the form (5.7). In 
this case we have three integrals, integration areas of which are defined by the appropriate (5- functions. Therefore 
we shall make transformations analogous to ( |5.9D . The transformations of the integration area of the second integral 
standing in (5.11) will be the same as (|5.9|). The transformation of the integration area of the third integral into the 



corresponding integration area of the first integral in (5.11) has the form 



h 



c 



(5.12) 



Ci = c 



,1 



/ / k * 

~~ k" ' 



C2 - C^77- 



Applying all of these transforms yields factorization of the integrand for the stationary equation ( 5.11 ) 



— / \Vk kl k 2 \ 2 &k-k 1 -k 2 5u-u 1 -n 2 Tk\k 1 k 2 dk\dk2 




Hence it follows that the figure bracket vanishes at a = — 5/2, (3 = — 2, namely, the solution has the form 

n k = Ak- 5 ^ 2 kl 2 , N k = Bk- b ' 2 kl 2 . (5.13) 

The obtained solution corresponds to spectra with constant energy flux P e , Connection between constants A and B, 
as before, is found from the stationary equation ( |5.l| ). From this equation it is possible to get the previous estimation 
on their ratio: c s A ~ VaB. 

The found above solutions of the Kolmogorov type are related to the only channel of interaction, namely, the inter- 
action between Alfvenic and slow magneto- acoustic waves that demanishes significantly the value of these solutions. 
Remind, that processes together with fast magneto-acoustic waves have growth rates of the same order of magni- 
tudes, and therefore they can not be ignored. Fortunately, the channel (with the fast magneto-acoustic waves) can 
be incorporated in the considered above scheme without essential generalizations. As was pointed out in the previous 
section, the maximal scattering of A-waves is attained at the maximal value of z-projection of momentum transfered 
to slow magneto-acoustic waves while scattering of A-waves. Therefore it is natural to assume that such behavior of 
the scattering amplitude of A-waves should lead to strongly anisotropic distribution of waves concentrated in a narrow 
cone of angles along mean magnetic field: k z 3> k±. Under this assumption it is possible to consider the dispersion 
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law of fast magneto-acoustic waves to be approximated by those for Alfvenic waves: u>2 ~ |& z |Va. Another important 
circumstance in this case is that the matrix element of interaction is diagonal with respect to polarization A: 



yAAi 

k'kt 



#AAi V k ,u k . 



Thus, for almost longitudinal (along the mean magnetic field) distribution there is almost no difference between 
Alfvenic and fast magneto-acoustic waves. Moreover, the direct energy exchange is absent between these waves. This 
means that for this region of angles Kolmogorov spectra for fast magneto- acoustic waves will have the same form as 
those of the obtained spectra ( 5.1C ) and ( 5.15 ). In this case in the expressions ( 5.10 ) and ( 5.13 ) Nk and B should be 



changed into Ni and B\, and the coefficient 



.4 



a-l/2 E g A 

1 EBx' 



The spectra, obtained in this section, will have the physical meaning if the locality property will be fulfilled. This 
requirement of locality consists in that contributions into interaction of the waves from both the pumping region 
and the dissipation region have to be small. The latter leads to the requirement of convergence of integrals in the 
equations ( |5.5| ) and fl5.6| ). 

Convergence of integrals relative to k z provides by the presence of two (5-functions containing k z . As far as con- 
vergence against transverse wave vectors concerns, the integrals are logarithmically divergent. The logarithmical 
divergence, to our opinion, is not so serious as a possible powerful one. Appearance of divergence is connected with 
bi-homogeneity of the probability U. If a medium would be isotropic and matrix elements V would would have 
the same degrees of homogeneity as for the MHD waves at [3 -C 1 (such situation, for instance, takes place for 
Mandelstamm-Brillouin scattering in isotropic dielectrics) then in such a case the locality property would be valid 
(compare with |34|). Violation of bi-homogeneity for interaction of the Alfvenic and slow magneto-acoustic waves 
appears for almost transverse propagation: k±/k z ~ /J -1 / 2 and for interaction of the fast magneto-acoustic waves for 
small angles ~ /3 1 / 2 . By this reason cut-off of integrals in the kinetic equations should be performed on the smaller 
angles: ~ /3 1 / 2 . Another possibility to avoid the logarithmic divergence is in seeking for solutions containing powers 
of logarithm from k± in ( 5.10| ) and ( |5.13 ). However, this procedure does not lead to determination of the powers, 
but, however, provides a convergence of the integrals. 

Last remark. The spectra ( |5.10 ) as well as ( [5.13 ) have the same power dependence on transverse momenta: 
nk,Nk ~ kj 2 . Their homogeneity degree against transverse momenta is the same as for 2D (5-function of k^. This 
means that, besides the anisotropic Kolmogorov spectra ( pM0| ) and ( |5.13 ), the singular Kolmogorov spectra are 
possible: 



n k = Ak z 2 S(k ± ), N k = Bk; 2 6(k ± ) 



and 



Ak~ b / 2 8{\<L A 



N k = Bk; 5 / 2 6(h ± ). 



Which spectra are realized indeed? Rigorous answer to this question is possible to get by stability investigation of 
the spectra or numerical experiment (in the latter case one has a hope on the qualitative understanding). Both these 
approaches require separate consideration. 
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